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Abstract 

This paper is about adiabatic transport in quantum pumps. The notion of "energy shift", a 
self-adjoint operator dual to the Wigner time delay, plays a role in our approach: It determines 
the current, the dissipation, the noise and the entropy currents in quantum pumps. We discuss 
the geometric and topological content of adiabatic transport and show that the mechanism of 
Thouless and Niu for quantized transport via Chern numbers cannot be realized in quantum 
pumps where Chern numbers necessarily vanish. 
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1 Introduction 

An adiabatic quantum pump |33j is a time-dependent scatterer connected to several leads. Fig. ^ 
is an example with two leads. Each lead may have several channels. The total number of channels, 
in all leads, will be denoted by n. Each channel is represented by a semi-infinite, one dimensional, 
(single mode) ideal wire. We assume that the particles propagating in the channels are non- 
interacting 1 and all have dispersion e(fc). For the sake of concreteness we shall take quadratic 
dispersion, e(fc) = k 2 /2, but most of our results carry over to more general dispersions. We also 
assume that the incoming particles are described by a density matrix p common to all channels: 
p{E) = (1 + e^- 6- ^)" 1 , with chemical potential p and temperature T. The scatterer is adiabatic 
when its characteristic frequency lo <C 1/t, with r a the typical dwell time in the scatterer. We 
take units so that ks = h = m = e = +1. 

An incoming particle sees a quasi-static scatterer. The scattering can therefore be computed, 
to leading order, by time-independent quantum mechanics, using the scattering Hamiltonian in 
effect at the time that the particle reaches the scatterer. In other words, we pretend that the 
Hamiltonian always was, and always will be, the Hamiltonian seen by the particle at the time of 
passage. This gives the "frozen S-matrix" . Since time-independent systems conserve energy, so 
does the frozen S'-matrix. We denote by S(E,t) the on-shell frozen S-matrix. 

The outgoing states have, to order lo°, the same occupation density as the incoming states 
(since S(E,t) is a unitary nx n matrix and the incoming densities are the same for all channels). 
This implies no net transport. However, at order lo, there is an interesting interference effect: An 
incoming particle of well-defined energy does not have a well-defined time of passage. This spread 
in time is a consequence of the uncertainty principle, and is not related to the dwell time of the 
particle in the scatterer. Thus, even if lo <C 1/t, the (frozen) S'-matrix seen by the tail of the 
wave packet will differ slightly from that seen by the head of the wave packet. This differential 
scattering causes the outgoing occupation densities to differ from the incoming densities to order 
lo , leading to a nonzero transport. 

The density matrix for the outgoing states, p ou t, is determined by p(Hq), the density matrix 
of the incoming states, and the S-matrix. Let Sd be the exact (dynamical) S-matrix. As we shall 
explain in section p ou t is given by 

Pout = p{H - £ d ), £ d ^iS d S2- (1.1) 

A dot denotes derivative with respect to time 2 . £d is the operator of energy shift introduced in 
|24j. It combines information on the state of the scatterer with its rate of change Sd- For time 

1 Quantum pumps have also been discussed in the context of Luttinger liquids, see e.g. |2ll30|. 
2 The time dependence of S<j is discussed in section 151 
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Figure 1: A model of a quantum scatterer with two leads. 

independent scattering, Ed = and Eq. p. 1(1 is an expression of conservation of energy. The 
formula for p ou t, Eq. p. 1(1 . holds independently of whether the scattering is adiabatic or not. 
We call the frozen analog of the operator Ed the matrix of energy shift. It is the nx n matrix 

£{E,t) =iS(E,t)S*(E,t), (1.2) 

and is a natural dual to the more familiar Wigner time delay |15j 

T(E,t) — —iS'(E,t)S*(E,t). (1.3) 

where prime denotes derivative with respect to E. As we shall see in section IY . 31 their commutator 



a = i[T,8] (1.4) 

has a geometric interpretation of a curvature in the time-energy plane, analogous to the adiabatic 
curvature. 

Adiabatic transport can be expressed in terms of the matrix of energy shift. For example, the 
BPT P3| formula for the expectation value of the current in the j-th channel, {Q)j takes the form: 

1 r°° 

(Qh(t)=-ir- / dEp'(E)£ jj (E,t). (1.5) 
27T Jo 

At T = 0, —p' is a delta function at the Fermi energy and the charge transport is determined by 
the energy shift at the Fermi energy alone. 

There are two noteworthy aspects of this formula. The first one is that (Q), which is of order 
w, can be accurately computed from frozen scattering data which is only an uj° approximation. 
The second one is that the formula holds all the way to T = 0, where the adiabatic energy scale u> 
is large compared to the energy scale T . 

The energy shift also determines certain transport properties that are of order uJ 2 . An example 
is dissipation at low temperatures. Let (E)j(t) be the expectation value of energy current in the 
j-th channel. Part of the energy is forever lost as the electrons are dumped into the reservoir. The 
part that can be recovered from the reservoir, by reclaiming the transported charge, is p(Q}j(t). 
We therefore define the dissipation in a quantum channel as the difference of the two. As we shall 
show in section 1531 the dissipation at T = is 3 : 

(E) 3 (t) - p(Q) 3 (t) = L^) ..^ t ) > 0. (1.6) 

Both the dissipation and the current admit transport formulas that are local in time at T = 0. 
Namely, the response at time t is determined by the energy shift at the same time. This is 
remarkable for at T = quantum correlations decay slowly in time and one may worry that 
transport at time t will retain memory about the scatterer at early times. This brings us to 
transport equations which admit a local description only at finite temperatures. 

The entropy and noise currents are defined as the difference between the outgoing (into the 
reservoirs) and incoming entropy (or noise) currents. Namely, 

Sj(t,fi,T) = s(p out ,j)-s(p j ), s(p) = — J dE(ho p){E,t) (1.7) 

3 For related results on dissipation see e.g. 1251 . For relations between dissipation and the S-matrix see 
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where [1311X2] 



J -rcloga; - (1 - a;)log(l - x) entropy, 

Kx) = \ f , s ■ ( L8 ) 

I x(l — x) noise. 



In the adiabatic limit, u> — > 0, and for w T \A t V r we find (see section l5~4"j) 

s i (*,/x,T) = /-Af|(M ! t)>0 ) entr ° Py; (1.9) 

27rfc J 16 noise, 

where 

A£] = (£% - (£ n ) 2 . (1.10) 

When T < uj the entropy and noise currents at time t are mindful of the scattering data for earlier 
times and there are no transport equations that are local in time. What sets them apart from the 
current and the dissipation is the non-linear dependence on the density. The non-linearity makes 
the transport sensitive to the slow decay of correlations. 

Our result about the noise overlap with results that follow from the "full counting statistics" 
of Levitov et. al. When there is overlap, the results agree. However, the results are mostly 
complementary, a reflection of the fact that both the questions and the methods are different. 
"Counting statistics" determine transport in a pump cycle in terms of the entire history of the 
pump. We give information that is local in time. In this sense, we give stronger results. On the 
other hand, the counting statistics determine all moments all the way down to zero temperature, 
while our results go down to T = for the current and dissipation but not for the entropy and noise. 
A detailed comparison of our results with results that follow from the Lesovik-Levitov formalism 
|21 | is made in appendix lAl For other results on noise in pumps see e.g., |27j . 

Transport in adiabatic scattering is conveniently described using semi-classical methods [22] 
a.k.a. pseudo-differential (Weyl) calculus [2E]- As we shall explain in section 03 S(E,t) is the 
principal symbol of the exact S-matrix, Sd- Semi-classical methods can be used to derive Eqs. (11.51 
11.61 ll.7fl . For an alternate point of view using coherent states see [S]. 

In section we discuss the geometric and topological significance of our results. We shall see 
that charge transport can be formulated in terms of the curvature, or Chern character, of a natural 
line bundle. This is reminiscent of works of Thouless and Niu |H4j which identified quantized charge 
transport with Chern numbers EE and inspired the study of quantum pumps. Nevertheless the 
situation is different here, since the bundle is trivial and, besides, the integration manifold has 
a boundary. This does not preclude the possibility that charge is quantized for reasons other 
than being a Chern number. In fact, one can geometrically characterize a class of periodic pump 
operations 7 for which the transported charge in a cycle, and not just its expectation value, is a 
non-random integer. It is to be cautioned that a small change of the scattering matrix will typically 
destroy this quantization. 



2 Pedestrian derivation of BPT 

At T = the Fermi energy is a step function and p'(E) = —S(fi — E), and BPT, Eq. (|1.5(l . takes 
the form 

(Q) j ( f M,t) = ^£ jj ( t i,t). (2.1) 

In this section we shall describe an argument [H] that explains this equation in the two channel case. 
The two channel case is special in that the changes in the scattering matrix break into elementary 
processes so that for each one BPT follows either from simple physical arguments or from known 
facts 4 . 

4 We assume that the transported charge depends only on S(fj.) and, linearly, on dS(fj,), regardless of the physical 
realization of the scatterer. 
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2.1 The two channels case 

In the two channel case, Fig. ^ the frozen, on shell, S-matrix takes the form 

S0>)=(l ( 2 - 2 ) 

with r and t (respectively r' and t') the reflection and transmission coefficients from the left (right). 
Eq. (0)1 reads 

2tt (dQ)_ = i(rdr + t'dt'), 2tt (dQ)+ = i(f'dr' + tdt) (2.3) 

(c?<5)_ is made from data (r, t') describing the scattering to the left, and similarly {dQ) + from 
those to the right 5 . 

To identify the physical interpretation of the differentials we introduce new coordinates (9, a, <f>, 7). 
The most general unitary 2x2 matrix can be expressed in the form: 

n (e ta cos0 ie-^sm6\ 

\ie^ sin 6 e" m cos6y { > 

where < a, <fi < 2ir, < 7 < ir and < 9 < ir/2. In terms of these parameters, the BPT formula 
reads 

2tt (dQ}± = ± (cos 2 9) da T (sin 2 6) d(f> - d<y. (2.5) 

As we shall now explain, the variations da, dtp and dj can be identified with simple physical 
processes. 

2.2 The snowplow 



Figure 2: Moving the scatterer changes a — > a + 2Uf d£ 

Let kp be the Fermi momentum associated with fi. Translating the scatterer a distance dt; 
multiplies r, (r') by e 2lkpd ^, (e~ 2 *' CjF ' £i ^) , and leaves t and t' unchanged. It follows that da — 2kpd^ 
corresponds to shifting the scatterer. 

As the scatterer moves, it attempts to push the kpdt;/ir = da/2ir electrons that occupy the 
region of size <i£ out of the way, much as a snowplow attempts to clear a path on a winter day. 
Of these, a fraction |t| 2 = sin 2 9 will pass through the scatterer (or rather, the scatterer will pass 
through them), while the remaining fraction |r| = cos 2 9 will be propelled forward, resulting in 
net charge transport of 

2tt (dQ)± = ± (cos 2 9) da, (2.6) 

in accordance with Eq. H2.5(l . 

Remark 2.1 It is instructive to examine the special case of a uniformly moving scatterer where 
we can use Galilei transformations to compute the charge transport exactly. By taking the limit of 
slowly moving scatterer we get a check on the result above. 

Since the mass of the electron is one, the Galilean shift from the lab frame to the frame of the 
scatterer shifts each momentum by — £. In the lab frame, the incoming states are filled up to the 
Fermi momentum kp while in the moving frame the incoming states of the ± channels are filled 
up to momenta kp ± £. In the moving frame, the outgoing states on the ± channels are filled up 

5 This is why S* is on the right in the energy shift, Eq. 11.21 . 
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to kp =p £, and partially filled with density |t'(fc)| 2 for momenta in the interval (hp — £, hp + £). 
Transforming back to the lab frame we find for Sp of the — (=left) channel 

'0 ifk<k F -2i 
Sp_(k 2 ) = { -| r '(fc + £)| 2 ifk F -2i<k<k F (2.7) 
if k > k F . 



Sp_{E) = -2k F i \T'{k F )\ 2 5{E - pi). (2.8) 
(Q)j(t) = ^- f°° dE S Pj (E,t), (2.9) 



2n Jo 



To order £, 
Since the current is 

Eq. \2. 6)) is reproduced. 

Remark 2.2 The net outflow of charge, (dQ)- + (dQ)+, vanishes to order £ but not to order £ 2 . 
This is because the moving snowplow leaves a region of reduced density in its wake. 

2.3 The battery 




A(x,t) 



o 



Figure 3: Applying a vector potential changes 4> ~ > 4> + J A 

To vary <j> we add a vector potential A. This induces a phase shift d<fi — J A across the scatterer, 
and multiplies t, (t') by e ldc ^, (e~ ld *), while leaving r and r' unchanged. This phase shift depends 
only on J A, and is independent of the placement or form of the vector potential. The variation in 
the vector potential induces an EMF of strength — J A — —<j>. To first order, the current is simply 
the voltage times the Landauer conductance |t| 2 /27r That is, 

2tt (dQ)± = T (sin 2 0) d<f>. (2.10) 

in agreement with Eq. I|2.5|l . 

Remark 2.3 Consider the special case of a time independent voltage drop. In a gauge where the 
battery is represented by a scalar potential, the pump is represented by a time independent scattering 
problem where the potential has slightly different asymptotes at ±oo. If the battery is placed to the 
left of the scatterer, the states of particles incident from that side are occupied up to energy p — <fi, 
while those incident from the right are occupied up to energy p. Suppose <f> is negative. Then 5p+ 
is 

5p+(E) = l\t{E)\ 2 ifp<E<p-4> 1 (2.11) 
[o tfE>p~<j> 
If, however, the battery is placed to the right of the scatterer then 5p+ is 

ifE < p, 

■ i ;;t(E+©)i 2 tfp<E<p~<j>, (2.i2) 

ifE>p-j) 
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In either case, to leading order in <fi, 



S P+ (E) = -<p\t(p)\ 2 S(E-p). (2.13) 
Plugging in Eq. \2.ty) . we recover Eq. fg.,51) . 

Remark 2.4 To leading order in <j), 5p+ is independent of whether the battery is to the left of 
the pump or to the right of it. To order (j) 2 this is no longer true as one sees from Eq. 
\2.Mjj) . The frozen S-matrix is, however, insensitive to the location of the battery. It follows that 
it is impossible to have a formula for Sp, accurate to order 0{uj 2 ), that involves only the frozen 
S-matrix and its derivatives (see ]2b^ for some model calculations in the non-adiabatic regime). 

2.4 The sink 

The scattering matrix depends on a choice of fiducial points: The choice of an origin for the two 
channels. Moving the two fiducial points out a distance d£ may be interpreted as forfeiting part 
of the channels in favor of the scatterer. This new, bigger, scatterer is shown schematically in Fig. 
0] This transforms the scattering matrix according to 

S(k F ) -> e ldl S(k F ), d 7 = 2k F d£. (2.14) 

This operation removes kpd^/ir electrons from each channel and so we get 

2vr (dQ)± = -2k F d£ = -dry, (2.15) 

in accordance with Eq. (|2.5|l . Changing 7 is therefore equivalent to having the pump swallow 
particles from the reservoirs. 




Figure 4: A scatterer that has gobbled up d£ of each wire. 

For arbitrary variations dS the above result still holds for the sum dQ_ + dQ + . This follows 
directly from a fact in scattering theory known as Birman-Krein formula |36j and in physics as 
Friedel sum rule ^B] which says that the excess number of states below energy p associated with 
the scatterer is (2ni)~ 1 logdet S(p), whence 

-2ir((dQ)_ + (dQ) + )(fi) = -idlogdetS(p) = 2d 1 . (2.16) 
2.5 The ineffective variable 

We have already seen that changes in a, ^,7 yield transported charges dQ± which are correctly 
reproduced by Eq. I|2.5|l . Moreover, for any change of s, the sum + (dQ) + is, too. To 

complete the derivation of Eq. H2.5H we must consider variations in 9, with a, <f> and 7 fixed, and 
show that (dQ)- ~(dQ) + =Q. 

Suppose a scatterer has a, 4> and 7 fixed, but 6 changes with time. By adding a (fixed!) vector 
potential and translating the system a (fixed!) distance 6 , we can assume that a = <j> = 0. Now 

6 This can be achieved, in general, only at a fixed energy, and we pick the energy to be the Fermi energy /x. 
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imagine a second scatterer that is the mirror image of the first (i.e., with right and left reversed) 
as in Fig. |3J Since 9 and 7 are invariant under right-left reflection, and since a and (j> are °dd 
under right-left reflection, the second scatterer has the same frozen S matrix as the first, and this 
equality persists for all time. 



N A 



Figure 5: An asymmetric scatterer and its image under reflection 

/From the frozen scattering data we therefore conclude that the currents for the second system 
are the same as for the first. However, by reflection symmetry, (dQ)- — (dQ) + for the first system 
equals (dQ) + — (dQ)- for the second. We conclude that (dQ) + — (dQ)- = for variations of 9. 

3 Alternative perspectives on BPT 

The pedestrian argument does not extend beyond the two channel case. This is because with 
more than two channels, a general variation dS cannot be described in terms of known physical 
processes. One can, nevertheless, derive BPT from general and simple physical considerations, 
without recourse to formal perturbation expansions in scattering theory. 

3.1 An axiomatic derivation 

The BPT formula for the current follows from the following natural axioms: 

1. Existence and Bilinearity: 

dQ m = o^kidSijSu + complex conjugate, (3-1) 

ijkl 

with universal (complex) coefficients a^ kl . 

2. Covariance: The formula is covariant under permutation of the channels. In particular, it is 
invariant under permutations of the channels other than the m-th. 

3. Gauge invariance: The formula is unchanged by time-independent gauge transformations, 
and also under time-independent changes in the fiducial points. 

4. Cluster: If the system consists of two subsystems, disconnected from one another, then 
the currents in each subsystem depend only on the part of the S matrix that governs that 
subsystem. 

5. Landauer: If a voltage, applied to a single lead vnl ^ to, is modeled by a time-dependent 
gauge transformation, then dQ m is given by the Landauer formula where the transmission 
probability is given by the scattering probability m' — > to. 

6. Birman-Krcin: 

YdQj = ^dlog(detS) = ±Tr (dStf). 

* — ' Z7T Z7T 

3 
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The physical motivations for most of the axioms are clear. For example, bi-linearity comes 
from the fact that the current is an interference effect between the original outgoing wavefunction 
(described by S) and an additional piece (described by dS). The one axiom that seems the most 
arbitrary is existence. Namely, the assumption that charge transport at T — is determined by 
the scattering matrix at the same time and at the Fermi energy alone. 

By covariance, it suffices to study the current on the first channel, dQi (and drop the superindex 
in a). Henceforth, Latin indices will run from 1 to n, while Greek indices will run from 2 to n. 

The most general bilinear (|3 - 1 1) vanishing identically on unitaries S and their variations dS is of 
the form a^u = ejiSik + dik5ji with hermitian matrices and (<ijfc). Indeed, by SS* — S*S = 1 
the matrices (dS)S* and (dS)* S are anti-hermitean, which implies vanishing of the stated bilinear 
form. We can thus fix ctijki in (|3.1|) by imposing the uniqueness constraints that ^ Q a a j a i and 
X)q a iaka are anti-hermitian. 

We shall now see that by axiom 3 

ciijki — CijSikSji (3.2) 
Indeed, if we move the fiducial point on the z-th channel by distance £j, the S'-matrix transforms 

by 

Sij - Sije^to+M. (3.3) 

and so Oijki — * fl^we**'*' . The only invariant terms are those with i = k and j = I or 
with i — I and j = k. Similarly, gauge transformations send 

Sij-tSye'te-M. (3.4) 

Now the invariant a%jki terms are those with i = k and j = I or with i = j and k = I. Putting the 
two selection rules together gives Eq. (|3.2|) . At this point Eq. (|3.1|) reduces to 

dQi — CijdSjjSjj + complex conjugate, (3-5) 

ij 

with the uniqueness constraints that cy are pure imaginary. 

By the clustering property c a p must all vanish, for otherwise dQi will be affected also by 
channels disconnected from it. Now c lQ = C2 and c al = C3 independent of a by permutation 
symmetry. We can therefore write Eq. I|3.I|) as 

dQi = cidSnSn + c 2 ^ <JSik§tk + c 3 ^ dSkiSki + complex conjugate, (3.6) 
k k 

Summing dQj over all channels, using that c\, c%, C3 are pure imaginary, and setting the result to 
agree with the Birman-Krein formula gives 

dQj = 2icilm ^2 dSjjSjj + 2(c 2 + c 3 )Tr (d55 f ) = ^-Tr (dSS^). (3.7) 

3 3 

Thus ci = and 47r(c2 + C3) = i. 

What remains is to distinguish between ci and C3 . Imagine modeling a voltage V on channel a 
by a time dependent vector potential that shifts the phase of the wavefunction on the a-th channel 
by dxj). Equating with Landauer gives 

dQi - -2z(c 2 |S lQ | 2 - c 3 |S Q i| 2 ) d<t> = 1- \S la \ 2 d<p. (3.8) 

Since (for n > 3) |5 a i| 2 and I^iqI 2 are independent, this implies that c 2 = i/Air and C3 = 0. We 
thus obtain BPT for n > 3 7 . 

7 To get the n = 2 case, consider a 2-channel scatterer as a degenerate 3-channcl scatterer, where the third channel 
is disconnected from the first two. Then 531 is identically zero, and the 3-channcl BPT formula for dQi reduces to 
the 2-channel formula. 
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3.2 Classical pumps 

The classical phase space associated to a given channel is the half-plane {x,p \x > 0, p € R}. We 
can choose coordinates so that points in phase space are labelled by the pair (E, t), where E is the 
energy of the (classical) particle and t its time of passage at the origin. For concreteness, let us 
assume a dispersion relation e(p) — e(— p) with e(p) increasing from to oo as p ranges over the 
same interval, e.g. a quadratic dispersion. Then 

E = e(p), t = -x/v, {v = e') (3.9) 

is a canonical transformation to the energy-time half plane {E,t \E > 0, t <G R}, since dE A dt = 
dx A dp. The mapping is singular when v — 0. States with t > arc incoming (at time 0), while 
those with t < are outgoing (at time 0). (All states are, of course, incoming in the distant past 
and outgoing in the distant future.) 

The phase space of n disconnected channels is T — U™ =1 I\ = {(E,t,i) \ E > 0, t 6 R, i = 
l,...n}. When analyzing pumps, i.e., channels communicating through some pump proper, T 
still serves as phase space of the scattering states. More precisely, (E, t, i) shall be the label of 
the scattering state whose past asymptote is the free trajectory with these initial data. Similarly, 
we may indicate a scattering state by its future oriented data (E',t',j). In this way we avoid 
introducing the full phase space for the connected pump. However, some of these trajectories may 
admit only one of the two labels, as they arc free for, say, t — + — oo but trapped as t — > +oo. With 
this exception made, the relation defines a bijective map, the (dynamical) scattering map: 

5:r-^T+, (E,t,i)»(E',t',j), (3.10) 

where T\r~ are the incoming labeled trajectories which are trapped in the future, and correspond- 
ingly for r \ r + . If (|3.10|) is viewed as a function of (E, t) with i fixed, the channel j is piecewise 
constant and the map to (E',t') symplectic. We shall illustrate S by an example in Sect. 16.41 The 
inverse map may be written as 

S- 1 : (E\ t', j) i * (E, t, i) = {E 1 - £ d (E', t',j),t' - T d (E', t', j), i), (3.11) 

which defines the classical energy shift Ed and the time delay Td as functions of the outgoing data. 
We remark that for a given static pump £ = and T is independent of t' . For adiabatic pumps, 
we have £ d (E',t',j) = 0{uj) and T d (E',t',j) = T(E',t',j) + 0{lo), where T is the time delay of 
the static scatterer in effect at the time of passage t 1 , on which it then depends parametrically. 
Since S^ 1 is volume preserving by Liouvillc's theorem, and the derivative w.r.t. time brings in a 
factor u>, we have 

£'+f = 0, (3.12) 

where £ is the part of £ d of order u . This relation shows that the static scattering data determine 
£(E',t',j) up to an additive function of t' , j and, as we shall see, cannot do better. This is in 
sharp contrast to the quantum case, where £ is fully determined by the frozen scattering matrix, 
see Eq. Ijl.2ll . We will further comment on the origin of this ambiguity in the classical case in 
Sects. I4~T1 16.51 and relate it to the lack of phase information in classical scattering. 

Similar to the quantum case, Eq. H1.5JI . is however the expression of the current in terms of the 
energy shift: 

/•OO 

Qj(t)=- dEg'(E)£(E,t,j), (3.13) 
Jo 

where g(E) is the phase space particle density in the incoming flow. We remark that in a semi- 
classical context g is related to the occupation density p by g(E) — p{E)/2-K. 

In fact, the net outgoing charge transmitted in the time interval [0,T] through channel j is 

Q 3 = f dE'dt' X[ o,T](t')g(E) - [ dEdt X[a .T](t)g(E), (3.14) 
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where E in the first integral is given through the map l|3.11[l if (E' , t'.j) e T+; if [E' , t',j) e T\r+, 
which may occur if E' is close to threshold energy and E < 0, we assume that g(E) = g(0), i.e., 
that the occupation of the bound states and threshold energies are equal. 

Eq. !|1.5fl is obtained immediately by expanding g(E) = g(E') — g'(E')£(E',t',j) + 0(u> 2 ) in 
the first integral 1|3.14[) . The contribution of the first term cancels against the second integral. 

Another derivation, which is more involved, is of some interest especially in view of the semi- 
classical discussion of pumps in Sect. 13.31 The first integral l|3.14|l equals 

r°° t=T 
dE'dt'x %T] {t)g{E) - / dE'g(E)T(E',t',j) . (3.15) 
r 3 Jo *=o 

In the adiabatic regime we may describe 1^ , to lowest approximation, as Tj = UjLjTy, where F^- 
consists of states (E',t',j) originating from lead i under static scattering. W.r.t. them and to 
next approximation, their preimages (E,t,i) appearing as arguments in the first integral (|3.15|) . 
are displaced by the vector field — (£(E',t',j),T(E',t',j)), which is typically discontinuous across 
the boundaries of the I\j's, but divergence free otherwise by (|3.12l) . (For an illustration, see 
Example 16.41 and Fig. 1101 there.) As a result, that integral differs from the second integral (|3.14|l 

by 



i=l 



(da E £(E',t',j)+da t T(E',t',j))xio,T](t)g(E), (3.16) 



where (d<jE,d<j t ) is the outward normal to STy. Within U^Lj^Fy we may distinguish between 
boundary parts contained in the boundary {E — 0} of Tj, and inner boundaries. The contribution 

of the former is dtg(0)£ (0, t, j) and, mostly for comparison with the promised semiclassical 
derivation, we formally write that of the latter as J r fl(E,t, j)X[o,T]{t)9(E), where fl(E,t,j) is a 
distribution supported on the inner boundaries. Since the map (|3.10() is bijective on T except for 
bound states, the displacements of the sets are such that YJj=i M(E, t,j) = 0. In summary, we 
obtain 



T 



t=T 




-+ 


f dE f 


t=0 


Jo Jo 



T 



Qi = J o dtg(0)£(0,t,j)- J dEg(E)T(E,t,j) _+ / dE j dtg(E)Q(E,t,j) (3.17) 
and, by differentiating w.r.t. T, 

Qj = g(0)£(0,t,j) - I dEg(E)T(E,t,j)+ I dEg{E)n{E 7 t,j). (3.18) 



The first term on the r.h.s. describes the release and trapping from bound states. The middle 
term describes the depletion of the outgoing flow as a result of a time delay increasing over time, 
since effectively no charge is exiting during a time dT. The last term describes electrons that are 
reshuffled between leads, but with no withholdings since Y^=i ^(E, t,j) = 0. 

^From (|3.18(l . Eq. (|3.13(l can be recovered: Let [E k (t), Ek+i(t)}, (k = 0, 1, . . .) be the intervals of 
the partition of [0, oo) on which £, T are continuous, and A/j£ = £(Ek + , t, j)— £{E^ — , t, j), Aj.T = 
T(Ek + ,t, j) — T(Ek — ,t, j), (k = 1, 2, ...) the values of their discontinuities at the endpoints. Then 

/>oo 

/ dEg(E) Q(E, t, j) = V g(E k )(A k £ - E k A k T), 
Jo k>i 

dEg(E)t(E, t,j) = -Y, dEg(E)t(E, t, j) + £ g(E k )E k A k T. (3.19) 



k>0 JEk fc>l 



Using Eq. I|3.12|l and integration by parts, the first term on the r.h.s. of l|3.19|l can be written as 



f ^ dEg(E)£'(E,t,j) = -g(p)£(0,t,j) -^ 5 (S fe )A fe £ f ^ dEg'(E)£(E, t, j). 

By collecting terms, we recover Eq. (|3.13p . 
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3.3 Currents and the T-S uncertainty 

We present a semiclassical derivation of Eq. I|1.5f) . where we however take for granted the physical 
meaning of energy shift of Eq. (|1.2fl . which will be established in Sect. 14. ll Since the time delay 
is a matrix, see Eq. Q1.3[l . we should, as a preliminary, point out that it is its diagonal element 
Tjj(E, t) which has the meaning of the average time delay of a particle exiting channel j. In fact, 
consider an incoming wave packet J dke~ l ( kx+e ( k ^ (p(k) in channel i centered on the trajectory 
—x = e'(k)t + c. The part of it scattered into lead j is J dke l ( kx ~ € ( k ^ Sji(e(k),t)<p(k) and is 
associated with x = e'(k)(t — (arg Sji)') + c (with ' = d/dE on argS^i), which implies a time delay 
of (arg Sji)'. Averaging this with the probability \Sji\ 2 for the particle to have come from channel 
i, we find for the average delay 

n n 

£ |S ri | 2 (argS>0' = Im^SjiS'^ = T n . (3.20) 

i=l i=l 

The net outgoing charge in the time interval [0, T] is at order u 

-i i>oo -i p oo /*T 

(Q)j = 7T dE ' dt'p(E)-— dE dtp(E), (3.21) 
27r Jo Jo 27r Jo Jo 

where E in the first integral is given through the map 

$ : (£', f ) ~ (E, t) = {E' - E i3 (E\ t'),t> - Tjj(E', t')) , (3.22) 

which describes the effect of the pump on the energy and the time of passage of an electron in 
terms of the outgoing data (E',t'). This is similar to the classical Eq. I|3.11|) except that £ and 
T are now defined in terms of the quantum mechanical frozen scattering matrix S, see Eqs. (|1.2I 
ll.3|) . Energies E' close to the threshold E' = may not be in the domain of the map $. Similarly, 
energies E may there fail to be in its range. These situations correspond to electrons released from, 
resp. trapped into, a bound state of the pump. 

The BPT formula, Eq. 1)1.5(1 . is again obtained immediately by expanding p(E) — p(E') — 
p' (E')£jj(E' ,t') + 0(uj 2 ) in the first integral (|3.21)l . The contribution of the first term cancels 
against the second integral. A further derivation, which is longer but adds another interpretation 
to the result, is by using $ as a change of variables. The Jacobian of l|3.22|l is 1 — fljj(E', i'), where 
Qjj is the divergence of the displacement (£jj,Tjj). As a matrix, f2 is the time delay-energy shift 
uncertainty introduced in Eq. (|1.4[> : 

n = i[T, £} = i (SS*' - S'S*) =£' + f . (3.23) 

Since is formally of order O(w), the Jacobian is close to 1 and the map (|3.22)l is invertible. After 
changing variables to (E, t) the first integral (|3.21|l extends over (E,t) £ $([0,oo) x [0,T]), which 
differs from [0, oo) x [0, T] to leading order by displacements —£jj(E, t) along E = and —Tjj(E, t) 
along t = 0, T. This yields 

r T oo r°° t=T r°° r T 

2tt(Q)j=- dtp(E)£jj(E,t) - dEp{E)T 3J (E,t) + dE dtp{E)VL 33 {E,t) , 
Jo Jo 4=0 Jo Jo 

and the differential version thereof is 

/•OO /'OO 

2%{Q) j (i) = p{0)£ 3J (0, t) - / dE p{E)tjj (E, t) + dE p(E) % (E, t) . (3.24) 



Upon inserting l|3.23|l and performing an integration by parts on /?£'„ , Eq. I|1.5|) is recovered. The 
interpretation of the three terms is the same as given after Eq. (|3.18|l . For the first term, with 
bound states now quantized, this is further discussed in the remark below. As for the last term, 
notice that X^=i ^jj = s ^ m holds true because of Eq. (|1.4)l . While generally, and in contrast 
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to the classical case, Q may have full support in phase space, it remains true that it vanishes if 
scattering is deterministic: If on some open subset of phase space 

S Ji (E,t)=0 for (3.25) 

where it is some fixed permutation of the channels, then S* is in the same manner related to 7r _1 , 
and £ , T to the identity permutation, i.e., they are both diagonal matrices. Hence f2 = by 1)1. 4[l . 



Remark 3.1 The first term on the r.h.s. of \3.24\l typically consists of delta functions located at 
times t where the pump has a semi-bound state at a band edge, i.e., a state which can be turned either 
into a bound state or a scattering state by an arbitrarily small change of the pump configuration. 
We illustrate this for e(k) — k 2 , and first claim: either S(0,t) = Iwheio S(E,t) = —1, or the 
(static) pump at time t admits a semi-bound state. This is seen as follows: For each k G C the 2n 
plane waves 

i smkx rooa-s 

eicoskx , e% — - — , (3. 2o) 

rZ 

({ e i}i=i being the standard basis of C n ), form a basis of solutions with energy k 2 in the n discon- 
nected leads. Upon connecting them to the scatterer an n-dimensional subspace of solutions is left, 
which depends analytically on k. Since the dependence of ]3.26\) is also analytic, solutions may be 
written as 

/ \ ^ ^ 7 1 Sill rvX r > r _\ 

Wk{x) = 2_^ a i e i cos kx + Pi e i — £ — > (3.27) 

1=1 

with amplitudes a — (a±, . . . , a n ), (3 = (fix, . . . , (3 n ) satisfying a set of linear equations 

A(k)a + B(k)(3 = (3.28) 



with analytic n x n coefficient matrices A, B. As i) defines an n-dimensional subspace, we 
have rank {A, B) = n. For k — 1^3. 2b}) reduce to 64, e^x. By a semi-bound state we mean more 
precisely a bounded solution ipo(x), i.e., one with f3 = in \3.2 r i\ . Its existence is tantamount 
to detA(O) = 0. If dctA(O) ^ 0, f^IW can he solved f° r a at smal1 k: a = ~-4(fc)~ 1 S(fc)/3, 
with arbitrary (3 G C™. For k > 0, the scattering matrix maps the incoming part of \3.2 r l\j to its 
outgoing part, ika + /3 = S(k 2 )(ika — [3). Thus 

S{k 2 )^{~ikA{k)- 1 B{k) + l){-ikA(k)- 1 B{k)~l)- 1 ~> -I , (fc -►()). (3.29) 

This proves the claim, and in particular that £(0,t) = except at times t when the scatterer has a 
semi-bound state. To discuss its behavior there, say this happens at t = 0, we assume generically 
that C{k, t) = B{k, t)~ l A(k, t) (= C(k, t)*) has a simple eigenvalue j(k, t) with a first order zero at 
k = t = 0. Let P be its eigenprojection at crossing and let a — sgn7(0, 0) be the crossing direction. 
Then, we claim, 

]im£(E, t) dt = -2naPS(t) dt , (3.30) 
so that in the process, by $2.1\) . the charge 

-Urn— / tv£(E,t)dt = a (3.31) 

is captured at arbitrarily small energy. Since the l.h.s. also equals lim^xo argdet S(E, t)|l £ , \3.tjl)) 
states that the phase of S(0,t) changes by 2tt upon capture of a bound state, which is a dynamical 
version of Levinson's theorem. The proof of $3.30)) rests on the approximation of , 



valid near k = t = 0. 



q(h 2 A _ ikP-W(0,0)k + j(0,0)t) 

b ^ k ifcP + (y( ,0)fc + 7(0,0)i) ' { ^ ' 
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4 Time dependent scattering and Weyl calculus 



4.1 The energy shift 

In this section we describe the notion of energy shift in the context of time dependent scattering 
theory and derive an operator identity relating the outgoing density p out to the incoming density 
and the energy shift. 

Energy is conserved in time independent scattering but not in time dependent scattering. This 
leads to an important differences between the S-matrix in the time independent and the time 
dependent case: In the time dependent case the S-matrix acquires a dependence on time shifts. 
Since this observation is not particular to pumps, it is simpler to describe it in general terms. 

Using the conventional notation of scattering theory, let Hq be our reference, time independent, 
Hamiltonian (often called the free Hamiltonian) and H = Ho + V the scattering Hamiltonian. V 
may be time dependent. We assume that H and H admit a good scattering theory. In the time 
independent case, conservation of energy is expressed as the statement that the scattering matrix 
S commutes with H (not HI). Hence, for the frozen S-matrix 

S fe - iHot = e- lHot S f . (4.1) 

This may be interpreted as the statement that the state ip, and the time shifted state e~ iHot ip 
both see the same scatterer. Therefore it makes no difference if the time shift takes place before 
or after the scattering. 

This is, of course, not true in the time-dependent case. The states ip and its time shift e~ tHot ip 
do not see the same scatterer. Consequently, the exact (dynamical) scattering matrix, Sd, for a 
time dependent problem, acquires a time dependence: 

S d (t)e- iHot = e~ iHot S d . (4.2) 

Now Sd{t) will, in general, not coincide with Sd except, of course, at t = since energy is not 
anymore conserved. In the time dependent case one does not have one scattering matrix, but rather 
a family of them, Sd(t). Since they are all related by conjugation, any one of them is equivalent 
to any other. To pick one, is to pick a reference point on the time axis. We shall use the notation 

s d (o) = s d . 

If Sd is unitary 8 , which we henceforth assume, then Sd = iSdS^ is Hermitian. We call it the 
"energy shift" for the following reason: S d satisfies the equation of motion 

iS d (t) = [Ho,S d {t)}. (4.3) 

Using the (assumed) unitarity of Sd, this can be reorganized as 

S d (t)HoS* d (t) = Ho-£d(t). (4.4) 

Conjugation by the scattering matrix takes outgoing observables to incoming observables. Eq. 14.4fl 
justifies identifying £ d with the operator of energy shift. 

Remark 4.1 If we let Qj denote the projection on the states in the j-th channel, then Qj = 
SdQjSd projects on the out states fed by the j-th channel. The energy shift generates the evolution 
ofQf. 

iQ = [£d,Qj] (4.5) 

We are now ready to derive Eq. ll.lj l which is an operator identity for p ou t- This is our 
starting point in analyzing adiabatic transport. By the functional calculus we can extend Eq. (|4.4|) . 
evaluated at t = 0, to (measurable) functions of Ho, namely 

Pout = S dP (H Q )S*d = p{H - £ d ) . (4.6) 

So far, no approximation has been made. The identity does not assume an adiabatic time depen- 
dence. 

^Sd is unitary as a map between the spaces of in and out states which may differ because states may get trapped 
or released from the pump. 
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Remark 4.2 For comparison we establish the classical counterpart to Eq. \4-4\l - Let T be the 
classical phase space described in Sect. VJ-'A and let (f> s : T — > T be the flow 4> s (E,t) = (E,t — s) 
generated by the Hamiltonian h(E,t) = E, i.e., the solution of the canonical equations of motion 

±Mx)=I(dh)\ Mxy (4.7) 

where x = (E, t) and I : T*T — > TV is the symplectic 2-form. If times of passage are measured not 
w.r.t. time but w.r.t. time t, then the scattering map, Sd(t), satisfies, cf. Eq. 

S d (t) °^ t =4>t° S d , (4.8) 

where 5^(0) = Sd has been introduced in Eq. $3.10\) . Since Sd(t) is a family of symplectic maps its 
vector field is Hamiltonian: 

pd(t)(x)\ t=0 = I(d£d)\ Sd(xy (4.9) 

where £d is a function on T uniquely determined up to an additive constant. By taking derivatives 
of I4.ti[ ) at t — we obtain I(d£d) + Sd*(I{dh)) = I(dh). Since Sd is symplectic we have Sd*I = I 
and hence Sd*(I{dh)) — (Sd*I)(Sd*dh) = I(S2)~ 1 dh = Id(h o S^ 1 ), so that we conclude 

hoS^=h-£ d , (4.10) 
provided the constant not determined by Sd is properly adjusted. 



4.2 The Weyl calculus 

A convenient language for discussing the relation between operators in quantum mechanics and 
functions on phase space, called symbols, is the Weyl calculus [221 128j . For pumps the classical 
phase space has been introduced at the beginning of Sect. 13.21 with points labelled by the pair 
(E,t), where E is the energy of the (classical) particle and t its time of passage at the origin. An 
additional index j = 1, . . .n labels the channels. 

The relation of a (matrix valued) symbol a{E, t) to the corresponding operator A is 



(t,j\A\t',f) = ±-J dEe~^ E a n , (e, 



t + t' 



(4.11) 



where \t,j) is the (improper) state in the j-th channel whose time of passage at the scatterer (in 
the Hq dynamics) is t. Equivalently, 

(E,j\A\E',f) =±-J dte^- E '» a jf (^,*) , (4.12) 

where \E, j) is the (improper) state in the j-the channel with energy E. It follows that (if A is 
trace class |3"T]1 

1 



Tr A = — / tradEdt (4.13) 

27T J 

where tr a denotes a trace over channels, i.e. a trace of finite dimensional matrices. Similarly, if a 
or b are (locally supported) functions, we have 

Tr (AB) = -!- / tr (ab)dEdt. (4.14) 

2-7T J 



5 Adiabatic transport 

The notion of approximation in adiabatic scattering requires some explanation. In this regime the 
scattering of a particle occurs on the time scale of the dwell time r which is short compared to 
the adiabatic time scale a; -1 . Therefore, the (unitary) operator Sd should be related to the frozen 
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scattering matrices S(E,t). While the uncertainty relation forbids specifying both coordinates E 
and t of a particle, the variables on which S(E,t) actually depends are E and Lot. This gives 
adiabatic scattering a semiclassical flavor where u> plays the role of h. Its theory can be phrased 
in terms of the Weyl calculus |221 128| with symbols which are power 9 series in to. In particular, as 
we shall explain, S(E, t) may be interpreted as the principal symbol of S d . The chain of argument 
in making the identification goes through the frozen S-matrix, Sf(t), where the time of freezing, t, 
is picked by the incoming state. 



5.1 Adiabatic scattering 

We shall show the following correspondence between operators and symbols which, in our case, are 
n x n matrix functions of E and t: 

S d (t ) S(E,t + to) + 0(u) (5.1) 

£ d <f=> £(E,t)+0(u 2 ) (5.2) 

Pout p(E)-p'(E)(S(E,t) + 0(LJ 2 )) + ~p"(E)£ 2 (E,t) + 0(w 3 ) (5.3) 

Since the Fermi function at T = is a step function, p' is a delta function and consequently, the 
notion of smallness in the expansion in Eq. (|5.3[1 is in the sense of distributions 10 . 

To see the first relation, let \t,j) denote the state that traverses the scatterer at time t and 
Sf(s) denote the frozen S-matrix associated with the Hamiltonian in effect at time s. Then 

(t,j\S d \t',f) = (t,j S f f,j') + OM (5.4) 

The matrix elements on both sides are significant provided t — t' is small, within the order of the 
dwell time, or the Wigner time delay. Using 



one finds 



(27T)- 1 / 2 J dEe lEt \E,j) (5.5) 



(t,j\S f (s)\t',j') = ^JdEe- i ^ E S(E,s), 8 = ^-, (5.6) 

and we have used the fact that Sf is energy conserving. This establishes Eq. (|5.1|l for to = by 
comparison with Eq. (|4.11(l . In the language of pseudo-differential operators S(E, t) is the principal 
symbol of S d . More generally, S(E,t + to) is the principal symbol of S d (to). The "quantization" 
of S(E,t) then satisfies Eq. Q4.2|) . as it must. 

Eqs. <|5.21 15. 3| ) now follow from the rules of pseudo-differential calculus [2H|, and the operator 
identity for the outgoing states Eq. (|1.1|) . 



5.2 Currents 

A rigorous derivation of BPT is presented in OEHl- Here, instead, we shall be content with a 
formal, but relatively straightforward derivation using Weyl calculus. 

Let Q l ^ out {x, tg) be the observable associated with counting the incoming/outgoing particles 
in a box that lies to the right of a point x in the j-th channel at the point in time to- The point x 
is chosen far from the scatterer, but not so far that the time delay relative to the pump is of order 
lo~ x . Namely, vt <C X <C v/lu. The symbol of Q m / out is a matrix valued step function: 

Q™ /out {xM)^ P j 6(v(t -t)-x)6(T(t -t)), (P j ) ik = 8j k 8i j , v = e'(p), (5.7) 

and Pj is the projection matrix on the j-th channel. Indeed, the position of a particle with 
coordinates (E,t) at time will be —v(t — to) at time to, see Eq. (13.91) . The particle will then be 

9 The dimensionless expansion parameter is lot. 
10 In agreement with Eqs. l2T8ll27l3l . 
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outgoing if t — to < 0. Notice that t = to falls outside of the support of the first Heaviside function. 
The associated incoming/outgoing current operators are 

Q i ; /ou \x, to) = i[H, Q" l/ou \x, to)] = i[H , Qj'° u \x, to)]. (5.8) 

Here we used the fact that beyond x, deep inside the channel, H(t) coincides with Ho- The symbol 
associated with the current is most easily computed recalling that in Weyl calculus commutators 
are replaced by Poisson brackets. This reproduces the usual notion of a current 

Q™ /out {x,to) ^ Pj {E, 0{v(t -to)-x)} = TP 3 S(t -to- x/v)e( T (to - t)), (5.9) 

x is where the "ammeter" is localized. By the assumption that the ammeter is not too far it leads 
to a slight modification of to, the time when current is measured. We henceforth drop x. Now the 
expectation value of the current is 

(Q)j(x,t ) = Tr {poutQT 1 ) + Tr ( P Q) n ) = Tr (Sp Q? ut (x,t )), Sp = p out - p. (5.10) 

Using Eq. I|4.14|l to evaluate the trace we find 

(Q)j(x,to) = / dEdtp'(E)S jj (E,t)5(t-t ) + O(u 2 ) 

= ~ j ' dEp'{E)E 3] (E,to) + 0(u 2 ) (5.11) 

reproducing Eq. I|1.5fl . 
5.3 Dissipation 

To compute the dissipation we start as in the previous section. Let £)" l / out denote the observable 
associated with the incoming/outgoing excess energy in the j-th channel in a box to the right of 
the point x. The excess energy is, of course, energy measured relative to the Fermi energy: 

D; l/out (x,t ) = \{Q l ; /out (x,tohHo-p} (5.12) 

and Q™l° ut ^ as m Eq. Ij5.7|l . The observable associated with dissipation current in the j-th 
channel is the time derivative of the excess energy, i.e., 

^.in/out/ , \ ^in/out/ , \ A.in/out/ , \ ■ jt ^.in/out/ , \ /rr io\ 

Dj {x,t ) = E j / (x,t ) - pQj [x,to)=t ffc-Dj (x,ta) ■ (5.13) 

The symbol corresponding to the dissipation current is then 

D) n/out {x, to) TPj (E - p) 5(t -t - x/v)0(T(t O - t)). (5.14) 

The expectation value of the dissipation current is therefore 

(D)j(xM) = Tr(5pD° ut (x,to)). (5.15) 

We shall now show that for T < \fw~fr the dissipation is quadratic in ui and is determined by 

Eq. fTjEJl . 

As in section l5?2*l we shall evaluate the trace using Eq. 1)4.14(1 . At low temperature p 1 is con- 
centrated near the Fermi energy. We may then approximate the energy shift up to its linear 
variation near p. For the term proportional to p' in the expansion Eq. (|5.3(l the contribution to 
the dissipation is proportional to 



1 

2^ 



dE P '(E) (£ n (p,t Q ) + (E- p)£' 3J (p,to) + 0{uo 2 )) (E - p) 

= 0(f3e-^)+0{uT 2 )+0{uj 2 T). (5.16) 
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The term proportional to p" in the expansion gives, 

i- J dEp"(E)(£ 2 ) n (E,t Q )(E-p) = ^(£ 2 ) rj {pM)+O{L0 2 T), (5.17) 

which is the requisite result, Eq. 11.6fl . 

The result 1)1. 6|) is remarkable in that we obtain the dissipation to order ui 2 by making two 
approximations, each valid only to order u). First we replace £d with £, and then we evaluate £ 
at E = p. Had we used this procedure to compute the quantities Ej and Qj separately, each of 
them would be off by nonzero 0(uj 2 ) terms (as can be seen in the snowplow and battery examples); 
nevertheless, the combination Ej — pQj is computed correctly to order u> 2 . The reason is that in 
each quantum channel one has the following lower bound on the dissipation [7|: 

Ej - pQj > ttQ 2 . (5.18) 

This bound is saturated by the outgoing population distribution that is filled up to energy p and 
empty thereafter. The dissipation should therefore be quadratic in the deviation of the outgoing 
distribution from this minimizer. Since the outgoing distribution is an 0{lo) perturbation of the 
minimizer, knowing the distribution to order ui should give the dissipation to order u> 2 . 

Remark 5.1 There would appear to be two problems with the argument above. In minimizing a 
functional on a region with a boundary, one obtains a quadratic estimate for the functional around 
its minimizer if the minimum occurs at an interior point. If the minimum occurs at the boundary, 
then a variation away from the boundary can increase the function to first order. Furthermore, 
whether the minimum occurs at an interior point or on the boundary, quadratic estimates depend 
on the Hessian being a bounded operator. If the Hessian is unbounded, then an arbitrarily small 
change in the point can cause an arbitrarily large increase in the functional. In our case, the 
minimum occurs at a point that is on the boundary of the constraints < p(E) < 1. There are 
large modes for the Hessian, involving adding electrons at arbitrarily high energies. 

Fortunately, neither exception is relevant. In fact, the correction of the distribution, as given 
by \5.$j) . consists of a local reshuffling of electrons around the Fermi energy and does not involve 
the large modes of the Hessian. These variations should not be viewed as being either towards or 
away from the boundary, since neither the expression 14 .b\ ) nor its opposite (replacing £ with —£) 
violates the constraints. 



5.4 Entropy and noise currents 

Entropy and noise introduce a new element in that the transport equation, Eq. I|1.7fl . depends on 
the density through a non-linear function h(p). For the entropy and noise h(x) is given in Eq. (|1.8|) . 
Using the fact that in either case h(0) = h(l) = 0, we shall show that for u> <C T <C yfujfr the 
currents are quadratic in uj and are given by 

Sj (t,p,T) = ^-A£ 2 (p,t))J^ dxh(x), (5.19) 

where A£ 2 has been defined in Eq. (|1.10|) . For the entropy the integral gives 1/2 and for the 
noise it gives 1/6. To complete the derivation of the noise and entropy currents, we now derive 

Eq. jnnnji. 

The condition oj <C T makes it possible to consider the outgoing state of the electrons with a 
fixed time of passage, provided the time resolution is short compared to uj^ 1 but large w.r.t. 7 7-1 . 
The state p ou tj — PjPoutPj is then given, see (|5.3|l . as 

PoutAE) = P(E - £jj(E,t)) + l -p"{E)({£ 2 ) 00 (E,t) - £ 33 {E,t) 2 ). (5.20) 

The entropy /noise current (|1.7|) is 

i- J dE((hop)(E-£jj(E,t))-(hop)(E))+±- J dE(h'op)(E)p"(E)A£ 2 (E,t). (5.21) 
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In these integrals, Sjj may be regarded as constant in E because of the condition T <C yoj/r. 
The first integral then vanishes and in the second we may pull AS J (/x, t) out of the integral. This 
leaves us with the integral 

f dE(ti o p)(E)p"{E) = -f3 f dph'(p)(l-2p) = 2/3 [ dph(p), (5.22) 
J Jo Jo 

where, in the second step, we have used a property of the Fermi function, p' = —(3p(\ — p), and 
an integration by parts in the last step. This establishes the result. 

We have nothing to say about the range T < ui. The noise at T = can be calculated using a 
formalism of Lesovik and Levitov that we discuss in the appendix. This formula is nonlocal in time. 
It is instructive to examine what goes wrong with our approach at T = 0. In this limit p' and p" 
are distributions and p a discontinuous function. Since it is not allowed to multiply distributions 
by distributions, or even by discontinuous functions, equations such as Eq. (|5.22|) which are not 
linear in p make no sense. This is a reflection of the fact that in the regime where T < u> entropy 
and noise currents have a memory that goes back in times of order (3. In the regime we consider 
to <C T the memory is short compared with the time scale of the pump and instantaneous formulas 
make sense. At the opposite regime, where T < u, a local formula in time cannot be expected. 



6 Examples 

Quantum pumps may be viewed either as particle pumps or as wave pumps. The particle interpre- 
tation has a classical flavor where the driving mechanism is identified with forces on the particles. 
The wave interpretation stresses the role played by phases and suggests that interference phenom- 
ena play a role. This duality can be seen in the BPT formula in the two channel case of section 
On the one hand Eq. (I2.5|l makes it clear that the phases in the S'-matrix, (a, <j), 7), play a role 
in charge transport. In fact, changing the transmission and reflection probabilities while keeping 
the phases fixed cannot drive a current. At the same time, the rate of change of two of the three 
phases, <j> and a, also admit a classical interpretation as EMF and Galilean shift. The pedestrian 
derivation is a reflection of the fact that particle interpretation is more intuitive. Here we shall 
consider two examples where dual reasoning is insightful. 



6.1 The bicycle pump 

In a bicycle pump the action of the valves is synchronized with the motion of the piston. The 
analogous quantum pump has synchronized gates as shown in Fig. The particle interpretation 
of the pump is simple and intuitive. The wave (or BPT) point of view is more subtle. In particular 
as we shall see, in terms of the elementary processes described in section the pump operates by 
changing the phases 7 and a: Galilean shifts arise from the synchronized action of the gates. 

Let us choose a length scale so that kp — ^ and an energy scale so that p = 1. In these 
units, choose the length of the pump, L, to be an integer L — n, pick the valves thin, S <C 1, 
and impenetrable, i.e., of height M with MS ^> 1. Consider the potential, shown in Fig. El that 
depends on two parameters, a and 6, that vary on the boundary of the unit square [0, 1] x [0, 1]: 



V a dx) 





aM 
106 

(1 - a)M 




if x < 0,, 
if < x < 5, 
if S < x < L 
if L <x < L + S 
if x > L + S. 



(6.1) 



Since the quantum box was designed to accommodate n particles, the pumps transfers n particles 
in each cycle. Like in the bicycle pump, at all times, at least one of the valves is closed. This give 
the particle point of view 11 . 



For related results see e.g. [S^ 
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aMn 



bM 



[1 



(l-a)M 



Figure 6: On the left a typical configuration of the potential of the quantum bicycle pump. On 
the right a loop in the a-b parameter space. 



We now consider this pump from the perspective of the phases in BPT. At all stages of the 
cycle, the transmission coefficients at the Fermi energy are essentially zero, so, by unitarity, the 
reflection coefficients r and r' are unit complex numbers. The phases of the reflection amplitude, 
7±a of Ea. H2.4(l . must therefore change by ±27rn in each cycle of the pump. How does this happen? 
As we shall see in spite of the fact that the pump is operated by manipulating gate voltages, the 
interpretation in terms of the ^-matrix is in terms of an interplay between Galilean shifts da and 
the Birman-Krein term dj. 

At a = 0, b = I , we have r = — 1, since the piston imposes a Dirichlet condition at x = 0, and 
r' = — exp(2ikpL) = — I , since the valve on the right imposes a Dirichlet condition at x = L. As b 
is decreased, r' remains —1 (since the right valve is closed) and by Eq. (|2.4() da = dj. Meanwhile, 
the wave functions of incoming waves from the left penetrate deeper and deeper into the region 
< x < L, eventually accumulating n half- wavelengths. The phase of r increases by 2im, since 
the left barrier has been effectively shifted a distance L to the right and / dQ\ = —n. 

The path with 6 = has no effect on the scattering matrix, since L is an integral number of 
half-wavelengths, so a Dirichlet condition at x — is equivalent to a Dirichlet condition at x = L. 
The left valve closes and the right valve opens but no current flows. The remaining two legs of the 
path can be similarly analyzed. Increasing b with a — 1 decreases the phase of r' by 2-7rn, while 
decreasing a with b = 1 has no effect on r or r'. 

The fact that t = t' = throughout the process might seem strange. After all, how can you 
transport particles without transmission? However, this is exactly what happens with macroscopic 
pumps. Good bicycle pumps are typically transmissionlcss, while bad pumps have leaky valves. 

6.2 The U-turn pump 

The U-turn pump, shown in Fig. is a highly schematic version of the quantum Hall pump. There 
are two leads connected to a loop of circumference I. The loop is threaded by a slowly varying 
flux tube carrying a flux $. The particle satisfies the free Schrodinger equation on the edges of 
the graph and and satisfies an appropriate boundary condition at the vertices. The boundary 
conditions are such that at the Fermi energy all particles are forced to make a U-turn at the loop. 
Namely, all particles coming from the right circle the loop counter-clockwise and the exit on the 
right while all those coming from the left circle the loop clockwise and exit on the left. 




Figure 7: The graph associated with a model of the quantum Hall effect. 
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Let us first look at this pump as a wave pump. Since all particles make a U-turn, the trans- 
mission amplitude vanish and the reflection amplitudes are phases. A left mover on the loop 
accumulates a Bohm-Aharonov phase in addition to the phase due to the "optical length" of the 
path. This means that the S-matrix at the Fermi energy is 



S(/x,$) 



By BPT the charge transport is 







(dQ)±(/i)=± 







2vr 



(6.2) 



(6.3) 



One charge is pumped from left to right in a cycle of the pump as $ increases by 2tt, the unit of 
quantum flux. 

The scattering calculation, although easy, does not really explain how the pump operates: How 
does it transport charges from right to left if all charges are forced to make a U-turn at the loop? 
The particle interpretation demystifies the pump: Particles in the loop see a force associated with 
the EMF $ which make the clockwise movers feel as if they are going uphill while the counter- 
clockwise movers all go downhill. Because of this some of the slow counter-clockwise movers turn 
into clockwise movers and exit on the other side. Although insightful, the particle interpretation 
does not readily translate to a qualitative computation without invoking some wave aspects. 

The integer quantum Hall effect can be described by a scatterer with four leads (north, 
south, east and west) with a north-south voltage and an east-west current. However, if the north 
and south leads are connected by a wire, and if the resulting loop is threaded by a time-varying 
magnetic flux to generate the north-south voltage, then one obtains a geometry shown in Fig. |H1 




Figure 8: The Hall effect as a pump driven by the EMF $. 



The U-turn pump of Fig. [7] models the essential features of this geometry. The 2D electron gas 
in the Hall sample, and the magnetic field applied to the Hall crystal, are modelled by the vertices, 
which scatter particles in a time-asymmetric manner. The clockwise movers of Fig. correspond 
to electrons that enter the Hall bar from the west, move along the edge of the crystal until they 
reach the south lead, go along the loop from south to north, move along the edge from north to 
west, and emerge to the west. The counter-clockwise movers correspond to electrons that go from 
east to north (along the Hall bar) to south (along the loop) and then to east and out the east lead. 
By standard arguments ^7], the edge states reflect the existence of localized bulk states in the 
crystal. 



6.3 A family of optimal pumps 

Optimal pumps [J] saturate the bound in Eq. (|5.18|) . In 0] it was shown that optimal pumps that 
do not break time-reversal are transmissionless. (The two examples above also fall into the category 
of being optimal and transmissionless.) The following example shows that a general optimal pump 
can have any value of r and t. 

In the battery of section l2~3"l the scatterer got in the way of the electrons, and the most efficient 
transport was with r = 0. In the snowplow of section l2~2"l the scatterer pushed the electrons and the 
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most efficient transport was with |r| = 1. In the following example of optimal pump we combine a 
voltage with a moving scatterer, such that the scatterer is moving along with the electrons, neither 
pushing nor getting in the way. In this case, the scatterer doesn't actually do anything, and we 
get efficient transport, regardless of the initial values of (r,t,r',t'). 

Write the scattering matrix of a system where a — 2/^£ and <fi evolve as 

S(l*) = te <* r / e 2 V? J , (6-4) 

Synchronizing the velocity £ with the voltage <p according to 

2(4 = <j) (6.5) 
makes S = iipa^S. The energy shift is then a diagonal matrix 

£ = iSS* = -j> a 3 (6.6) 
which implies that the pump is optimal. 



6.4 The phase space of a snowplow 

We give a description of a classical snowplow moving on the real axis at speed vq during the 
time interval [— T, T], but at rest before and after that. It is described by the (total) phase space 
R 2 3 (x,p) with Hamiltonian function h — p 2 /2 + V(t), where V(t) is a barrier of fixed height 
V > Vq/2 and zero width located at 

(-vqT, (t<-T), 
x(t) = lvot, (-T <t < T), (6.7) 

(v T, (t>-T). 

We use the notation of Sect. 13.21 and denote by 1 the left channel (x < 0) and by 2 the right one 
(x > 0). Let C Tj C T+ be the outgoing labeled trajectories originating from channel i and 
eventually ending in channel j. The same meaning has T~ C T, C T - , except that the trajectories 
are incoming labeled. 

If a particle crosses the above scatterer of zero width, then its trajectory is free at all times. 
Hence 

i7 3 = r+, r 2 - x =r+. (6.8) 

It suffices to compute these subsets only. The remaining ones are then given by complementarity: 

r i - 1 ura = r i , r+ur+^r,, (6.9) 

with disjoint unions. 



• r i2 = ^12 ■ Depending on their time of passage t, trajectories of this type will (see Fig. 
require an energy E = v 2 /2 > E c = w 2 /2, with critical energy E c given as 



E c = V if \t\ > T ■ 
resp. by v c - v a = \fW , i.e., 



f'o 



T = T 



vq r 
J2E~ 



1 < ^ -2 -f Ul s T V ° X. 



E c = -(V2V + v o y if \t\<T- 



This portion of phase space is drawn dark shaded in the upper part of Fig. 1101 
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Figure 9: The thick line represents the position x(t) of the snowplow. The other lines are free 
trajectories of common energy E indicated by the slope v — \J2E > 0. Their times of passage t 
are read off at the intercepts with the abscissa. 



• r 
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T^. In this case the slope of free trajectories is v 



critical energy E c is 



E r = V 



E c = ~{VW - v ) 2 



if 
if 



1*1 >T 



JTe 



/ 2E < 0. As a result the 

T, 
T. 



2E 



This portion of phase space is drawn dark shaded in the lower part of Fig. ^| 



6.5 Classical scattering from a battery 

This example shows that, in the classical case, static scattering data cannot determine the en- 
ergy shift. Consider the classical version of the battery, see Sect. 12.31 with Hamiltonian function 
h(x,p) = (p — A) 2 /2 and gauge A(x,t) — t(f>'(x) of compact support. Clearly, this describes parti- 
cles which get accelerated as they cross the pump, whence there is an energy shift. More formally, 
the quantity f(x,p) = (p — A) 2 /2 + (f> is a constant of motion, as it is verified from Newton's 
equation (d/dt)(p — A) = —A = —<fi'. In the leads, / = p 2 /2 + </>, which implies that the energy 
E = p 2 /2, as defined there, gets shifted by £ — — (fil 00 ^ if a particle crosses the battery from left 
to right. On the other hand, for each static scatterer, p — A is a constant of motion, whence the 
static scattering maps all equal the identity map! The quantum mechanical phase information, 
which was present in the static ^-matrix and determined the energy shift, is of course unavailable 
here. 



7 Geometry and topology 

When a pump goes through a cycle, so does the scattering matrix S(E,t). The rows of matrix 
define unit vectors in C" and the charge transport in the j-th channel can be interpreted in terms 
of geometric properties of these vectors. 

As we shall explain in section [7"2l the charge transport is naturally identified with the global 
angle accumulated by a row of the S-matrix during a cycle. It need not be an integer multiple of 
2ir. Computing this angle is formally the same as computing Berry's phase JOj- That the global 
phase has direct physical significance is related to the fact that a quantum pump is a wave pump. 

Interestingly, the charge transport in a closed cycle can also be computed by forgetting al- 
together about the (global) phase provided one has knowledge about what happens on a surface 
spanning this cycle. This reflects a basic geometric relation between the failure of parallel transport 
and the curvature. 
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Figure 10: Portions of phase space corresponding to transmitted (dark shaded) and reflected (light 
shaded) trajectories with labels +/— corresponding to outgoing/incoming data. The curves on the 
right halves correspond to t = T ± -2=8= T. 



7.1 Charge transport and Berry's phase 

Let 

(Sji\ 



\1>i 



s 



•j2 



G C" 



(7.1) 



\SjnJ 

be the transpose of the j-th row of the frozen S matrix, evaluated at the Fermi energy. Yet another 
rewriting of the BPT formula at zero temperature is as 



2ir(dQ)j — i(ipj\dipj 



(7.2) 



The expression on the right hand side is familiar from the context of adiabatic connections and 
Berry's phase |10j . To simplify notation we shall fix a row, say the first, and drop the index j. 

There are important conceptual differences between the Berry's phase associated to a quantum 
state and the phases that arise in the study of a row of the S-matrix \ip). In the usual Berry's 
phase setting one starts with a circle of Hamiltonians to which one associates a (unique) circle of 
projections l^)^, say on the ground state. To represent these projections in terms of a circle of 
eigenvectors one needs to choose a reference phase arbitrarily for each point on the circle. The 
physical adiabatic evolution picks its own phase relative to the reference phase. Berry's phase then 
measures the phase accumulated in a cycle 12 . In particular, it does not make sense to talk about 
the phase accumulated on a path that is not closed. 



By choosing the energy of the state to be zero one can always get rid of the dynamic phase. 
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In contrast, a cycle of the pump is a circle of vectors of the scattering matrix (not a cycle 
of projections). No choice of a phase needs to be made — the phases are all fixed by solving the 
scattering problem. In particular, the vector \tp) always returns to itself after a complete cycle of 
the pump. Moreover, unlike the usual Berry's phase, Eq. l|7.2[l makes perfectly good physical and 
mathematical sense also for an open path and not just for a closed cycle. 

What, then, does the Berry's phase measure in the present context? This is addressed in the 
next section. 

7.2 Global angle 

As we now explain the charge transport in a quantum pump has a geometric interpretation in 
terms of the global angle of \ip). 

The notion of global angle is obvious in the case that the circle of vectors are all parallel, i.e., 

m = ^ |Vo> (7.3) 

with \ip ) fixed. The global angle is 7 and is related to the charge 2w(dQ) = —dj. 

Now, how do we compare global angles when vectors are not parallel? This is a classical 
problem in geometry whose answer relies on the notion of parallel transport. Namely, one needs a 
rule for taking one vector to another without changing its global phase. A natural way to do so is 
to impose that there be no motion in the "direction" \ip)(ip\. Explicitly, we say that is parallel 
transported if 

M<#ty)=0 (7.4) 

We can now compare the global phase for any two vectors connected by a path. The path may 
be either open or closed, and the global angle makes sense in either case. The right hand side of 
Eq. H7.2JI is (minus) the change in the global angle and the left hand side identifies it with the 
transported charge. This identifies charge transport with a global angle. 

The global angle need not change by a multiple of 2n in a cycle of the pump — except in the 
special case that — e %1 1 V'o) with |^q) fixed. The failure of parallel transport for closed paths 
is interpreted in geometry as curvature. 

7.3 Curvature 

In the previous section we identified charge transport with the global angle of \ipj)- Remarkably, 
one can compute the charge transport in a closed cycle while forgetting altogether about the global 
phase and relying only on the projection {ipj\. For a closed cycle one can relate the line integral 
on the boundary of a disk dD to a surface integral on the disk D via Stokes 

2<!r(Q) j =i[ (^dt/jj) =i [ (#j|<%> (7.5) 

JdD J D 

In the context of pumps the identity is known as Brouwer's formula 11 . 

Remark 7.1 For other ways to rewrite Eq. |7.5| ), use the identities 

i dipj) = -i(dS A dS*)jj = -iTrPjidPj A dP^Pj, (7.6) 

where Pj — S*PjS is the projection onto the state feeding channel j . The r.h.s. is the trace of the 
curvature of the connection Pjd (see e.g. *5J, Sect. 9.5). 

We shall now describe a different interpretation of Brouwer's formula that focuses on the Wigner 
time delay and the energy shift. We consider the charge transport in a cycle, so t is an angle. By 
Eq. H3.24f> . and assuming no semi-bound states, so 8(0, t) = 0, the charge transport in a cycle is 

2w{Q)j = / dE AdtQjj (7.7) 
Jc 
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where C is the cylinder [0, fi] x S. Now 

Qjj dE Adt = dSjj A dt + dE A dT j:j = -i (dS A dS*) r] . (7.8) 

The difference between this formula and Eq. (|7.5|l is the domain of integration: A disc in Brouwer's 
formula, and a cylinder here. However, since £(0, t) = the bottom of the cylinder may be pinched 
to a point and the cylinder turns to a disk. The r.h.s. is the trace of the curvature of the connection 
Pjd (see e.g. jS], Sect. 9.5) or of its connection 1-form i[(dS)S*]jj = Ejjdt — TjjdE. Similar 
equations are found in the context of the quantum Hall effect [32|, see e.g. [S] Sect. 11.3), but, 
unlike there, Jt § dEAdt Q(E, t)jj is not a Chern number as a rule, since the integration manifold 
has a boundary. 

7.4 The two channel case 

The two channel case is particularly simple. = (^rj lives in S 3 . The projection associated to 
\ip) can be identified with a point on S 2 according to 

W (VI - 1± tP (7-9) 

where h is a unit vector in R 3 and a is the triplet of Pauli matrices. Eq. I|7.5[) then says that the 
charge transport is half the spherical angle. 

The claim that one can compute the charge transport in a closed cycle while forfeiting all 
knowledge of the global phase means, in the 2-channel case, that all that matters is z — r/t'. z 
lives in CP 1 = C U {oo} and is related to n above via the usual stereographic projection 

|V.) - (2Rc(rF),2Im(rF),|rP- |tf ) = (f^,^,^) , (7.10) 

The curvature 

j(drdr + dt'dt') = i(l + \z\ 2 )- 2 dzdz, (7.11) 

can be written entirely in terms of z. The current pumped by a small loop in parameter space can 
therefore be computed by a calculation on S 2 . 

For large loops, however, things are more subtle. If a path in parameter space maps to the 
equator, is the bounding region the northern hemisphere (area +27r) or the southern hemisphere 
(signed area —2tt)1 Because of this ambiguity, the path in S 2 only determines the fractional part 
of the expected charge transported, not the integral part. In particular, if the path on the unit 
sphere is trivial — the ratio z — r/t' is constant — then the fractional part is zero and the charge 
transported is an integer. The bicycle pump, where t' is identically zero, is an example of this 
quantized transport. 

How, then, can one determine also the integral part without a knowledge of the global phase? 
Consider a region in parameter space bounded by our path. The integral part is how many times 
it wraps around S 2 (i.e., the degree of the map). For the bicycle pump, this is counting how many 
points on the unit square [0, 1] x [0, 1] correspond to a particular value of z, say z = Q. In fact, 
since we know, be elementary arguments, that the charge transport in each cycle is n, there must 
be (at least) n points in the interior of the square where the potential is reflectionless at the Fermi 
energy. 

7.5 Chern and winding numbers 

Part of the motivation for pumps is as standards of charge transport, whence quantization is an 
issue. In this section we want to explain why Chern numbers in quantum pumps must vanish. 
Instead, in some cases, winding numbers play a role. 
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The Hopf fibration gives the unit sphere S 2n 1 C C" the structure of a U(l) principal bundle 
over CP™" 1 via the map 

n . S 2n-l > CP" 1 (7.12) 

W — IMVI, (7-13) 

where we realize CP"" 1 as the set of all rank-1 projections in C™. 

Chern numbers in transport theory |35| typically arise as follows. The physical parameter space 
is a closed surface, A4, say a sphere or a torus. The Hamiltonian, acting on the vector space C™, 
is then a function on A4, and any one of its simple eigenvalues map M. to CP" -1 . The pullback of 
the Hopf fibration is then a U(l) bundle over parameter space, with local geometry and possibly 
nontrivial topology. The integrated curvature is a Chern number, and may be nonzero. 

In adiabatic pumps, however, the frozen S matrix defines a map from A4 directly to S" 2 ™ -1 , and 
indirectly to CP™ . This gives a trivialization of the bundle, and shows that all Chern numbers 
are zero. Besides, as seen from Eq. (J73J, charge transport is the integral of the Chern character 
i (dipj\ dipj) over a surface D with boundary, and thus not given by a Chern number. 

More appropriate is the point of view on pumps taken in the previous section, or its general- 
ization to n leads. If a pump operation is of the form 17.3fl . then §{dQ) is the winding number 
of e %1 |2I Ej . Alternatively, the curvature Eq. (|7.6[) can be computed in terms of Pj = Tr(\ip}), cf. 
Eq. (|7.11|) . and the integral l|7. 5fl thus performed over tt(D). The condition 17. 3|) means dir(D) — 0, 
which again shows that §{dQ) is an integer. 

An example where this quantization occurs is the U-turn pump, section T6. 21 The charge trans- 
port in a cycle is the winding number of the map from the circle of fluxes to the circle of complex 
numbers of modulus one. A small change of the parameters underlying Fig. |H1 will only modify 
in the effective description by Eq. (|6.2|) and thus preserve quantization. On the other hand, 
if the U-turn pump is viewed as an example of a Schrodinger operator on a graph, then a generic 
perturbation in this class will destroy quantization. 

7.6 Geometry of dissipation and noise 

Eq. H1.6J) and Eq. ljl.9|) that describe dissipation and entropy (and hence also noise) currents have 
a simple geometric interpretation in terms of the fiber bundle 5 2 ™" 1 — > CP™" 1 . The j-th row of 
£ describes the velocity of \ip) in S* 2 "" 1 . Of this, £jj is the projection of this velocity onto the 
fiber, and £jk, with k ^ j, give the projection of this velocity onto CP"" 1 . The current (Qj), 
and the minimal dissipation \£jj\ 2 /Air, are both functions of motion in the fiber, while the excess 
dissipation 

(E s ) - li{Qj) - *{Q 3 ) 2 = i- l^-l 2 (7-14) 

*" Mi 

is the "energy" (that is, squared velocity) associated with motion in the base. In particular, a 
pump operation is of the form l|7.3|l if and only if the bound Eq. I|5.18|l is saturated, or equivalently 
if the noise at T > or at T — 0, see Eqs. (|1.9IIA.19|l . vanishes. Such pumps may be called optimal 
w.r.t. the j-th channel. 

In an interesting piece of work Mirlin and Makhlin 23 , relate the problem of finding a cycle 
with minimal noise production at T — to the problem of finding a minimal surfaces supported 
by a given loop. This result is outside the scope of transport properties which are local in time, as 
it deals with T = and the noise. 
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A Comparison with the theory of full counting statistics 



The formulas for the entropy and noise currents, Eq. (|1.9I) . are singular in the limit T — ► 0. This 
appears to be in conflict with the Lesovik-Levitov's formula (LL) which gives finite noise at T = 
(see below). We shall verify here that, in fact, for T > u>, LL is consistent with our results. We 
shall start by recalling LL. 



A.l The Lesovik-Levitov formula 

Here we shall describe a version of the Lesovik-Levitov formalism where LL is an identity 
rather than an (adiabatic) approximation. 

Assume that in the distant past and distant future the Hamiltonian H(t) of the pump coincides 
with Hq, the disconnected pump. Let Qj be a projection on one the j-th channels. Since the channel 
is fixed we suppress the index j below. It is important in this formulation that the channels are 
disconnected at the distant past and distant future and that [Q, Ho] = 0. It is also assume that 
the initial state of the system is thermal state p{Hq). The counting statistics can be described 
by means of the characteristic function x(A) = J^^L-oc PnG lXn , where p n is the probability for n 
charges having been transferred to channel j in the course of whole process. The formula is 

X (\)=det(l + p(e- lXQ / 2 e^e- lX Q/ 2 -l)^, Q = S* d QS d , p = p(H ) (A.l) 

where, as before, p is the Fermi function and Q is a projection on the states feeding the channel 
in question. 

Remark A.l The reason Q plays a special role is traced to the fact that the second quantized 
dT(Q) is the number (charge) operator in the channel. 

Since Q is a projection e ^ = 1 + (e lA — 1)Q. Using this one finds 

= iA sin A - (1 - cos \)A 2 
+ 2isin(A/2)(l - cos(A/2))(QQQ ± + Q ± QQ) 



e -*AQ/ 2 el \Q e - t \Q/2 _ 1 = ^ a s in A — (1 — cos X)A 2 (A.2) 



where 

A = Q-Q (A.3) 
is a difference of projections. To second order in A 

logx(A) = Tr\og(l + i\pA-^pA 2 ^j + 0(A 3 ) (A.4) 

- i \Tr(pA)-^Tr(pA(l-p)A)+0(\ 3 ) 
In the last step we used the fact that Q commutes with Ho and the cyclicity of the trace. 

A.2 Charge transport 

The expectation value of charge transport into the j-th channel is the first moment: 

(Q)j = -i(logx)'(0) = Tr (pA) = Tr ((p(H - 8 d ) - p(H ))Qj) (A.5) 

This is clearly the correct result, independent of the adiabatic limit, for the rhs is precisely what 
one means by the change in the total charge in the j-th reservoir. In the adiabatic limit £ d is small 
and hence 

(Q)j « -Tr ( P '{H )E d Q ) =~J dtdE p\E)8 n {E,t) (A.6) 
in agreement with BPT. In the last step we used Eq. (j4.14|l and Eq. 1)5. 3JI . 
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A. 3 Splitting the noise 

Noise is the variance of the distribution associated to xW 0T i more precisely, the variance per 
unit time. It splits into two positive terms. One term is proportional to the temperature — this 
is the Johnson-Nyquist noise. The second term involve correlations at different times and survives 
at T = 0. This is the quantum shot noise. 
The variance is 

((AQ) 2 ) = -(log X )"U=o = Tr ( P A(1 p)A), (A.7) 

Now, write 

P A(l - p)A = P A 2 (l - P ) + P A[A, p] = P (l - P )A 2 + p[p, A] A. (A.8) 
Using the cyclicly of the trace and the average of the two terms in Eq. 1|A.8|) , we find 

Tr(pA(l-p)A) = Tr(p(l-p)A 2 )+±Tr([p,A][A,p}). (A.9) 
Each term is positive. The Johnson-Nyquist noise is the first term 

Qjn = Tr (p(l - p) A 2 ) = — TTr (p'A 2 ) > 0, (A.10) 
and the quantum shot noise is the second term: 

Qqs = \ Tr {[p, Q] [Q, p}) = \ Tr ([Sp, Q] [Q, Sp}) > (A.ll) 

(We have repeatedly used [Q, Ho] = 0). Since the semi-classical limit of a commutator is of order 
H, the quantum shot noise vanishes in the classical limit. 

A. 4 Thermal noise 

p' is a multiplication operator in E which, at low temperatures, is concentrated near p. The symbol 

Aj ^> a(E,t) = S*(E,t)PjS{E,t) - Pj. (A. 12) 



of Aj is 



where Pj is an n x n matrix that projects on the j-th channel. The Johnson-Nyquist noise at low 
temperatures can be written as 

2nQ 2 JN = t[ dtty(a 2 (p,t)) = T [ dttr (P(t) + P - P(t)P - PP(t)) 
Jr Jr 

= 2T]T / dt\S\%(p,t), (A.13) 

where tr denotes a trace of n x n matrices, and we have used the fact that P and P are one 
dimensional projections. We see that the Johnson-Nyquist noise is proportional to the temperature 
and the time integral of the conductance at the Fermi energy. It is finite since for large times H{t) 
coincides with H$ and the scattering matrix reduces to the identity. 

A. 5 Shot noise at finite temperatures 

The symbol associated to [Sp, Q] is, to leading order, the matrix 

- P '{E)[E{E,t),Pj\ 

where, as before, Pj is the projection on the j-th channel. At finite temperatures p'{E) is a smooth 
function and it makes sense to look at the square of the symbol as we must by Eq. (|A.11|) . Now 

tr [£(E, fyPjliP^Sfrt)] - 2((£ 2 ) n - (£ n ) 2 )(E,t) = 2A£^(E,t). (A.14) 
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Substituting in Eq. (|4.13|) and using the fact that p 1 is localized near the Fermi energy, we find 

= ^Jdttefat) (A.15) 
in agreement with Eq. I|1.9|l . 
A.6 Shot noise at T = 

The result at finite temperature may suggest that the noise over a pump cycle diverges as T — > 0. 
This is not the case. At T = the symbol p' is a distribution and it is not permissible to multiply 
them as we did in Eq. (|A.15(1 . However, the problem can be easily avoided by simply not using 
the r.h.s. of Eq. IjA.lljl but instead the middle identity. As before, we approximate the symbol of 
Q by its value on the Fermi energy 

Qi S*(E,t)PjS(E,t) « S*(n,t)PjS(n,t) = q 3 (t) (A.16) 

This approximation makes Qj a multiplication operator in t and hence 

Qqs= I I dtdt' p(t-t')\r(q(t)-q(t')) 2 , (A.17) 

JR JR 

where p is the Fourier transform of the Fermi function. At T — the Fermi function associated 
with chemical potential p is 

'<*> = ^Tho) 6 " 1 "- (A ' 18) 

The shot noise at T — 0, in the limit of large p, is given by 

Iff 1 - \ (S(p,t)S*(p,t')) ? 
QlsM = ±>} J dt dt ' ^ (A.19) 

Since S is unitary the numerator vanishes quadratically as t — t' — > 0. It follows that the integrand 
is a bounded function. However, it is supported on set of infinite area made of two strips: One 
along the t axis and one along the t' axis. Nevertheless, the decay properties of the denominator 
one easily sees that the integral is convergent. The noise in a pump cycle does not diverge as 
T ^ 0. 
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